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The continuous spontaneous localization (CSL) theory of dynamical wave function collapse is an 
experimentally testable alternative to non-relativistic quantum mechanics. In it, collapse occurs 
because particles interact with a classical random field. However, particles gain energy from this 
field, i.e., particle energy is not conserved. Recently, it has been shown how to construct a theory 
dubbed "completely quantized collapse" (CQC) which is predictively equivalent to CSL. In CQC, 
a quantized random field is introduced, and CSL's classical random field becomes its eigenvalue. 
In CQC, energy is conserved, which allows one to understand that energy is conserved in CSL, as 
the particle's energy gain is compensated by the random field's energy loss. Since the random field 
has energy, it should have gravitational consequences. For that, one needs to know the random 
1 field's energy density. In this paper, it is shown how to construct a symmetric, conserved, energy- 

momentum-stress-density tensor associated with the quantized random field, even though this field 
obeys no dynamical equation and has no Lagrangian. Then, three examples are given involving the 
random field's energy density. One considers interacting particles, the second treats a "cosmological" 
•^T ' particle creation model, the third involves the gravity of the random field. 

^ ■ PACS numbers: 03.65.Ta, 03.70.+k, 02.50.Ga 
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I. INTRODUCTION 



The CSL (Continuous Spontaneous Localization) theory[l| describes wave function collapse as a physical, dynamical, 
process. A motivation for thus altering standard quantum theory is that CSL has a well-defined interpretation. By 
comparison, no one has been able to graft a well-defined interpretation onto the formalism of standard quantum 
theory. A crucial test for such an interpretation is that it give rules for specifying the realizable (i.e., the possible) 
states of nature, and their probabilities of realization. There are no such rules for standard quantum theory: people 
decide upon the possible states in ad hoc fashion, utilizing educated judgments about specific physical situations. 
■ In CSL, the Schrodinger equation is modified by adding a term containing a classical white- noise field u>(x, i). 

Every state evolving under a possible field is a possible state. However, the theory also provides a rule for applying a 
• probability measure to each random field. The result is that only a subset of random fields evolve probable states. The 
dynamics is such that a superposition of states, which differ sufficiently in their mass density distributions, rapidly 
evolves to one such state. That is, the theory describes the macroscopic world we see around us. In many cases, 
CSL provides the same experimental predictions obtained from standard quantum theory, but there are experiments 
which can distinguish between themf^]. As emphasized by BohmQ, this possibility of suggesting and stimulating 
experiments, which might not otherwise be envisaged or undertaken, provides a strong reason for exploring the 
consequences of an alternative to quantum theory, were there no other reason. 

One experimentally testable consequence of CSL is that particles gain energy, e.g. electrons in atoms and nucleons 
in nucleii should "spontaneously" be excited. Experiments have show that the white-noise-particle coupling has to be 
mass-proportional, or otherwise a discrepancy would have been observed. This is suggestive of a connection between 
collapse and gravity, which a number of authors have proposed for other reasons 

Energy non-conservation in a theory is not considered to be a good thing. However, recently, a theory called 
CQC (Completely Quantized Collapse) has been investigated which, it has been shown, is completely equivalent to 
CSL[fj|. CQC has in it a quantized random field W(jx., t), which commutes with itself everywhere in space-time, and 
its eigenvalues are the CSL w(x,t). Since it is a Hamiltonian theory, energy is conserved: the increase in particle 
energy is precisely compensated by the decrease in the random field's energy, and this is then true of CSL as well. 

But, with these new ideas, there is the potential for new physics. A new source of energy ought to have new 
gravitational effects. However, to calculate those effects requires knowing the energy density. The purpose of this 
paper is to supply an expression for the energy density. 
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In Section II we define W(x, t). Section III defines its energy-momentum. Section IV (and Appendix A) presents 
a symmetric conserved energy-momentum-stress density tensor whose spatial integral equals its energy-momentum. 

Section V gives a brief introduction to CSL, with pertinent equations. Section VI contains a discussion of CQC, 
and explains how conclusions about CSL can be drawn from CQC, in particular, how to express the stress tensor in 
the language of CSL. 

Section VII and VIII discuss collapse dynamics for non-relativistic, interacting particles. Expressions for the 
ensemble-average particle and random field energy-densities are given. Section IX gives the comparable expressions 
for a model of'cosmogenesis, with a Hamiltonian which creates particles out of the vacuum. 

Section X presents an example where the random field energy density interacts gravitationally with the particle 
mass-density. Collapse dynamics is considered for a superposition of static particle mass-density distributions. The 
result is that collapse occurs as usual, except that th e state ve ctor at time t describes the state of the particles on 
the proper time hypersurface t' = t[l + </>(x')] ~ ty/— ffoo(x') (f, x' are the space-time coordinates, (f>(x') is the 
gravitational potential, and goo( x ') = — 1 — 20(x') is the 00th component of the metric tensor in the limit of a weak 
gravitational field). 



II. QUANTIZED WHITE NOISE FIELD 



We define a quantized white noise field 



W(x) ee ^ \ dk[e ik -*b(k) + e- ik -*tf(k)], (1) 



(2tt) s 

where A is a constant, £ M = (x, t), — (k, k°), k ■ x = k • x — k°t, and 

[&(*), &t(jfe')] = 6(k-k'), [b(k),b(k')}=0, [tf(k),tf(k')}=0. (2) 
We also define a conjugate field 

n(x) = 2 ^ J dk[-e ik -*b(k) + e-^bHk)]. (3) 

It easily follows that 

[W(x),Il(x')]=iS(x-x'), [W(x),W(x')] = 0, [Il(x),11(x')}=0. (4) 

For fixed x, W{x) has an eigenvector with eigenvalue w(x) which can take on any value between — oo and oo. 
Because W(x) has the unusual property for a quantized field that it commutes with itself everywhere in spacetime, a 
joint eigenvector of W(x) for every x can be constructed, 

W(x)\w) = w{x)\w), (5) 

Since w (x) can take on any value at any x, it is a sample classical white noise field. 

The vacuum state |0) satisfies J dkf(k)b(k)\0) = 0, where f(k) is an arbitrary function. Therefore, from |T]) and 
©, [W(x) + 2i\U(x)}\0) = 0, or 

(w\W(x) + 2iAn(ac)|0) = [w(x) + 2XS/5w{x)](w\0) = 0, (6) 

( 5/5w(x) is the functional derivative). The solution of (J6j) is 

{w\Q) = e- ( -* xrl ^-" dmo ^ x) . (7) 

The limits on integrals over x in this paper are always (— oo, oo), so an explicit limit, as in Eq. ([7]), always refers to t. 
The states (w\ are summed over utilizing Dw ~ \\ x dw(x), where the constant is chosen so that Dw\w)(w\ = 1 
(a rigorous definition can be given by discretizing space-timeQ), e.g., J Dw(0\w)(w\Q) — (0|0) = 1. 
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III. W-FIELD ENERGY-MOMENTUM 

The four- momentum, which generates space-time translations of W(x), II(:r), is 

PI ee J dkk v b\k)b{k). 

By solving Eqs.JI]), © for b(k), we find 



b(k) 



1 



(2^ 



dx 



1 



2VA 



w(x) + i-An(x) 



— ik-. 



(8) 



(9) 



Substituting 
there results: 



into ([8]), using expik ■ x = —i(d/dx^)expik ■ x in the expression for b>, and integrating over k, 
d 



P" 



—i / dx 



dx h 



1 W 2 {x) + ^IL 2 (x)} + f dxl\W(x)-^-U(x) - U(x)-^-W(a 
2 J ./ 2 ra„ - ' 



8A 



dx- [W(x) — U{x) - U(x) — W(x) 

da;W(a:) — H(ar) = - / (MI(x) — W(a;) = - dx —W(x) TL(x) 
ox v J ox v J vax v 

dx{W(x),-VW(x)}Il(x). 



(10a) 
(10b) 
(10c) 
(lOd) 



In going from (10a) to (10b), we have assumed no contribution from the boundary terms. This may be achieved 
e.g., by defining W(x), H(x) to vanish outside a suitably large space-time hypervolume. 

In (10c), we also recognize that the operators can be commuted freely within the integral, since dS(x — x)/dx u = 0. 
Eq.(lOd) displays {P°,P W }. 

IV. W-FIELD ENERGY-MOMENTUM-STRESS DENSITY TENSOR 

We want to construct a stress tensor T£"(x) for the VF-field, consistent with the constraints that it is symmetric, 
satisfies d v T£ u {x) — and P^ = f dxT° v (x). It is well known that stress tensors are not necessarily unique. For 
example, one can add a term to a field theory Lagrangian, of the form of a space or time derivative of a function 
of fields, which does not change the equations of motion but which does change the canonically constructed stress 
tensor. Therefore, we shall be content to display a stress tensor which satisfies these constraints. 

Appendix A contains a "derivation" of this stress tensor expression. Indeed, there appears an added term which 
is dropped without violating the above constraints. Not only does this dropped term make no contribution to the 
energy-momentum, but it creates and annihilates "W-particles," whereas the T£"(x) which remains and is displayed 
here does not. 

Without further ado: 



T^{x) 



— r= / dx\dx 2 [ dk\dk 2 5(k 2 — k 2 ) sQ(ski)Q(sk 2 ) sin[fci • (x — x\) — k 2 ■ (x — x 2 ) 
> n y j j - ^ 



(11a) 



2(2tt) 

■[k^k v 2 + k\k» + (1/2)77^(^1 - k 2 ) 2 ]W(x 1 )n{x 2 ) - n{x 1 )W(x 2 )} 

7 / dx\dx 2 I dk 1 dk 2 8(k 2 — k 2 ) sQ(ski)Q(sk 2 ) sin[fci • (x — X\) — k 2 ■ (x — x 2 )] 
\ J J s=-l 

7T-7T- + 7T-7T- " + /-) 2 [W( Xl )U(x 2 ) - U( Xl )W(x 2 )}, (lib) 

oxi/j, ax 2v dx 2fi ax\ v ox\ ax 2 

In (fTTjl . is the step function and the metric tensor's non-vanishing diagonal elements are rf-^ = (—1, 1, 1, 1). 

We note that T^ v (x) = T^ M (x). Now, consider d v T£ v {x). The derivative, acting on the sine, results in a factor 
[k\ v — k 2v \. Then, factors in the integrand are 



8{k{ - k 2 )[k lu - k 2y ][k^k v 2 + + (l/2)r 7 ^(fc 1 - fc 2 ) 2 ] = 5{k 2 - k 2 )(l/2)[k^ + U 2 l ][k 2 - k\ 
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so d v T^ v (x) = 0. It is remarkable that this equation describing the local flow of four- momentum is satisfied, given 
the very non-local nature of T£ v {x) displayed in Eqs. (|13Hi"5")) below. 
Finally, consider the spatial integral of (11a): 



J dxT°» = Jd Xl dx 2 y - & 1 dfc 2 ^5(A ; °-^)^se(<)e( S fc°)5(k 1 -k 2 ) s in[fc 1 -( 



2(27 

■2k° 1 k^[W(x 1 )n{x 2 ) - U{ Xl )W(x 2 )} 



x 2 - x\. 




cos[fci • (xi - x 2 )][W(xi)n{x 2 ) - Ii{xi)W{x 2 )) 
( d 



dx\ v 



IL( Xl ) \W{x l ) 



(12a) 
(12b) 

(12c) 



which is the result (10b), so = J dxT®"(x). Thus, all three constraints are satisfied. 
In obtaining (12a), the integral over x gives a factor <5(ki — k 2 ), so 



<J(ki - k a )e(sjfe?)e(«Jfe§)5(fc? 



kl) 



J(ki - k a )e(*A?)$(fc? - k Q 2 )/2\kl\ 



Thus ki — k 2 , so the term ~ rf v vanishes, and \k\k 2 

In going from (12a) to (12b), we particularly note use of the identity s9(sfc°)/|fc°| = 9(sfcJ)/fc°. 

In going from (12b) to (12c), we first integrate by parts, moving the derivative onto the operators, then integrate 
over k\ which results in 6{x\ — x 2 ), and then integrate over x 2 . 

The stress tensor expression (lib) contains the form factor 



1 2 J — ^ny-^ fv^ . 



G(x — xi, x — x 2 ) = —r= [ dkidk 2 5(k\ — k 2 ) s@(sk1)Q(sk 2 ) sin[fci • (x — x%) — k 2 ■ (x — x 2 )\. 
(2tt)' J ■t-f 



(13) 



The stress tensor is not relativistic, because of the s-dependent terms in (|13p . nor is it a local expression. However, 
(|13p is responsible for the relativistic and local nature of the four-momentum, since, following the steps outlined for 
Eqs.Cr 



J dx.G(x-x 1 ,x-x 2 ) = J dfc 1 dfc2^5(fc 1 -fc 2 °)^ S e( S A : °)e( S fc 2 ) ( 5(k 1 -k2)smhfc 1 .( a ; 1 -X2)] 



5(xx - x 2 )^y J dfc?-^ sin[fc?(ti - t a 



= j<y(xi - x 2 )e(*i - t 2 ) 



(14) 



i.e., the time derivative of (TT4")) is the local Lorentz scalar. 
G is calculated in Appendix B: 



G(x — x%, x — x 2 ) 



if 



(2tt) 4 da 2 



T x Q(a + if) T 2 e(- ( J + T 2 2 ) 



(15) 



where a ■ = (x — xi) 2 — (x — x 2 ) 2 , Tj = t — ti, and "P denotes taking the principal part. Its non-locality is evident in the 
location of the events xi, x 2 where a vanishes, whose neighborhoods make the largest contribution to G at the event 
x. If one constructs the forward and backward light-cones at x, and considers the family of (hyper-) hyperboloids of 
two sheets within these cones and asymptotically tangent to them, as well as the family of hyperboloids of one sheet 
outside these cones and tangent to them, a vanishes when both Xi and x 2 lie on the same hyperboloid. However, in 
the application we shall be considering, this non-locality is severely truncated by considerations of what is probable, 
as opposed to the above adumbration of what is possible. 



V. CSL 



The CSL Schrodinger-picture evolution of an initial state vector \<p) is defined as 



(16) 
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In (|16p . T is the time-ordering operation, H p is the particle hamiltonian, and 

- J^WI J d*e-^^ 2 ±^M n ti( Z )U*), (17) 

where the mass density operator (the sum in (|17j) ) is expressed in terms of the creation and annihilation operators 
at x of the nth type particle of mass M n , and M is the proton mass. A and a are respectively a collapse rate and a 
mesoscopic distance, parameters introduced by Ghirardi, Rimini and WeberQ for their collapse model, provisionally 
given the values A -1 « 10 16 sec, a w 10~ 5 cm (but, if the theory is proved correct, ultimately to be determined by 
experiment). 

The evolution (|16|) is non-unitary, so the state-vector norm is not preserved. CSL's rule for the probability density 
to be assigned to w(x, t) is 

p„=2(v,#,«)2 (is) 

with Dw = Jl x t dw(x, t)/ ' ^2irX/dxdt (rigorously definable in discrete space-time), so that / DwV w = 1. 

Eqs. (|16HT5| completely define CSL. With H p = 0, if states \aj) describe different mass density distributions, it can 
be shown that an initial state \ip,0) = J^j c j\ a j) evolves into ~ |aj) as t — > oo, provided 

lim T" 1 / dtw(x,t) -> 2A / (a,-|A(x)|a,-). 
Jo J 

Moreover, it follows from (fT5|) that the probability measure of all such states is \c } \ 2 . All fields with different asymptotic 
time averages have probability measure 0. With H p ^ 0, this same behavior obtains provided the collapse dynamics 
is rapid compared to the Hamiltonian dynamics, which is usually the case. 

Expressions simplify in the interaction picture, where the state vector only evolves due to its collapse dynamics. 
Defining the interaction picture state vector \ip,t) w = expiH p t\il), t)f,, (TH))) becomes 

\i/),t) w =Te-ix/o dt '/^'h(x'/)-2AA(x',t')] 3 |^_ (19) 

where A(x, t) = exp(iH p t)A(x) exp(— iH p t) (i.e., the particle creation and annihilation operators in (JTTJ) become time 
dependent operators). The probability expression <| 1 8|) is unchanged, except that the <S"s are removed. It follows from 
this revised (fP8f and (fl9|) that the density matrix which describes the ensemble of state vectors which evolve under 
all possible w(x, t) is 

p{t) = /'.Diu|i/>,f)t £ «,(V',*| 

= Te"^o *t'f<te'[AL(x\t')-A R ( x \t')]* fQ\ (20) 

(p(0) = \(f>)((/)\). The first expression in (20) is so simple because the normalization factors which divide the un- 
normalized state vectors cancel the probability factor {T5| . In (f2"0"|) , the operators with subscript L (R) appear to the 
left (right) of p(Q), and T time-orders the operators to the left, and reverse-time-orders the operators to the right. 



VI. CQC AND CSL 



The Schrodinger evolution of the state vector in CQC is 

|*,i) S = e -"[«P+ P °+2A/dxA(x)n(x,0)]|Q^|^_ / 21 s 

where P° is given by (lOd), II(x, 0) is given by (|3|), and (w\0) is given by 0. The CQC interaction picture state 
vector is = exp it[H + P°] |*, t) s , which yields 

(w\y, t) = Te~ l2X /o dt f dx^(x,t)(5/«t«(x,t)) e -^ JZe dt' J dx't« 2 (x',t') |^ 

= 7- e -Tx/-°°oo dt '/' ix '["'( x '.*')-2A J 4(x',t')e(t-t')e(t')] 2 |^ (22) 

using |J7]) ■ One way to think of this is to regard W(x, t) as a "pointer" located at x, labeled by t (so there 
is an infinitude of pointers at x). The wave function for this pointer has an initial (broad!) gaussian distribu- 
tion exp — (l/4:X)dtdxw 2 (x,t) of possible "positions" w(x,t). It evolves only during the time interval (t,t + dt) to 
exp—(l/4X)dtdx[w(x,i) — 2XA(x, t)] 2 , describing a standard (inaccurate!) local measurement of A(x,t). 
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The CQC state vector (|22|) is identical to CSL's (fT9]l . except for an extra numerical factor f(t): 

(tu|tf,t) = f{t)\i>,t) w where f{t) = e~™ f °°° dt ' /<**'™W) e -£ JT df /dxV(x',t'). 

/(i) does not play a role in calculations. It describes the wave functions for inessential pointers which cither never 
will make measurements (t' < 0) or which have yet to make measurements (t' > t). Following the usual quantum 
mechanical rule, the probability density associated with the eigenvalue w(x,t) of W(x) is | {w\^, t) | 2 which, when 
integrated over the inessential values of w(k' , t'), is identical to the CSL probability rule (fT8)) . So, CSL and CQC are 
completely equivalent when it comes to making predictions concerning particles. 

However, CQC does more. It naturally provides an expression for the energy P® of the white-noise field, as well 
as a particle-field interaction energy, such that the total energy H = H p + P^ + 2X J (fx'A(x')II(x', 0) is conserved. 
And, as we have now seen, it allows definition of a stress tensor associated to the white-noise field. 

The reason these structures have appeared in CQC, and not in CSL, can be seen by considering the CQC expectation 
value of any hermitian operator functional F{W(x),H(x)} which is polynomial in its arguments (but which may 
depend upon particle operators in any way), where the time argument of x lies in the range (0, i), and then expressing 
it in terms of the CSL state vector: 

(V,t\F{W(x),Il(x)}\y,t) = J Dw{V,t\w){w\F{W(x),IL(x)}\y,t) 

s c s 

Dw(*,t\w)F{w(x),—^-}](w\*,t)= Dw w (^,t\F{w(x),—--}\^,t) w . (23) 



iSw(x) ' J ' ' iSw(x) 

The functional derivative which appears in this expression means that F is not a quantity which can be attributed to 
a single CSL state vector \ij),t) w . F"s expectation value depends upon a set of CSL state vectors infinitesimally close 
to \tp,t) w . At first glance, this appears to conflict with the accepted ontology of CSL (which is perhaps why no one 
looked for such expressions). That is, reality is described by a single state vector. Nature chooses a single w(x) to 
evolve the state vector which is actually realized, and all other w(x)'s represent possible but unrealized evolutions. 

However, two considerations may stimulate allowing such expressions within this ontology. One is that, for state 
vectors with infinitesimally close w(a;)'s, any difference in their description of nature is indiscernible. The other is that 
CSL already has consideration of an ensemble of infinitesimally close w(x)'s fundamentally built into it. According to 
(I18p . Dw w (ip,t\ip,t) w is the probability to be associated with an ensemble of state vectors which have evolved under 
fields in an infinitesimal neighborhood of w(x). Thus, one may think of the family of state vectors infinitesimally 
close to the realized state vector as belonging to that state vector. 

Therefore, consider the CSL expression 

w ^,t\F{w(x), -^}\ij,t) w + w y,,t\F{w(x), ^£^}M,t) w 

(24) 



\,{ip,t\i),t) w 

(the arrows show the direction in which the functional derivative operates). This is real and, when it is multiplied by 
the probability Dw w (ip,t\^,t) w and w is integrated over, with use of integration by parts and hermiticity of F, the 
result is the ensemble expectation value (|23j) . According to the argument given here, (|24[) may be regarded as the 
expectation value of the operator F for the CSL state \tp,t) w . 

In the remainder of this paper, we shall discuss the ensemble-average energy-momentum density in CSL, utilizing 
the CQC formalism to do so. 

VII. ENSEMBLE-AVERAGE ENERGY-MOMENTUM DENSITY FOR PARTICLES 

Consider particle operators in the Heisenberg picture. The annihilation operator for a particle of mass m is 
£(x) = £(x, t) — exp(iH p t)^(x) exp(—iH p t). For free non-relativistic particles, where H p = (2m) -1 J dxV^ (a;)-V£(x), 
there is not an energy-momentum-stress density tensor. There is a mass density tensor T^ u (x), 

T°°(z) = m£l(x)Z(x), T%(x)=T™(x) = ±{-[d£(xMx)+liHx)ld i ax)}} 

2,1 

T%(x) = (Amr'i-id^ix)}^)} - [d^(x)m(x)] + mzHxMW+tHxMd&x)]} (25) 
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whose component (mass-density flux) is the momentum density, and there is an energy density tensor Tg u {x), 



T e 00 (x) = (2m)- 1 Ve t W-Ve(x), T°\x) = Tf{x) = ^{-[^ i^Wi^x) + d k d^ ( x )[d k Z(x)}} 

T%(x) = (2m)- 2 {-[d k ^(x)][d k d i d j ax)] - [d k d t d^(x)][d k ax)} 
+ [d i d k ^(x)][d j d^(x)} + [djd^ix^did^ix)}} 



(26) 



whose r°° component is the kinetic-energy density. The free particle tensors in Eqs. (|25l26p are conserved, 9 (1 T'"' = 0. 
[There is a conserved free particle relativistic energy-momentum density tensor, with components 

T ao (x) = 2- x { [V-V 2 + m^{x)] £(x) + ?(x) [^-V 2 + m 2 £(x)] }, T 0l {x) = T%(x), 



T lJ {x) 



^(x)d i d j ^(x) + [d'&^ixMx) - [^(x)]^^) - [d J ^(x)]d^(x) 



(V 2 acts on £'(x) and V| acts on ^(x)) which we shall not use here.] 

Suppose the particles interact through a potential, so H p includes the potential energy term 



V p = 2' 1 J d^ 1 d^(x 1 )^(x 2 )V(^ 1 - xa)t(xx)t{x 2 ) 



Since the potential is non-local, one cannot define a local energy density. One may define an energy density which 
depends on non-local contributions, 



H p {x)^T™(x) + 2-^\x) 



iff/ 



J dx^t^O^X^x)^!) 



which has the virtue that its spatial integral equals H p . However, one cannot form a symmetric conserved tensor with 
Ti v {x) as its 00th component. The mass-density tensor is also not conserved when there is an interaction: 



d u T^(x) = [l-S^(x) 



fa), 



i.e., the divergence is the force density. Nonetheless, we shall define the energy-momentum density operator for the 
particles as Pg(x) = (H p (x), l%(x)). 

Now, consider the particle energy-momentum density in CQC. Its expectation value (equal to the CSL ensemble 
average) for an initial particle state \<f>) and the initial vacuum field state |0) is, with help of (21), 

%{x) = s (*,t|^(x,0)|*,t) 5 

= (4>\ (0\e tt[H " +p ^ +2X f WWll^fx, 0)e~ tt[H " +p ° +2A f dxA ( x ) n ( x >°)]|0)|</!>) 

\Te i2X f° dx ' A (x'Wx')-p^ x ^ e -i2\ f* dx'A(x')H(x') |^ ^27) 



We note that (27) is expressed in terms of Heisenberg operators for the particles. We now continue, by inserting 
J Dw\w)(w\ into (27), and utilizing (22): 



V^x) = (4>\T / rHue-w^^^C^-aAA^pM^e-^/o^'Itu^-a^C*')]'' 



(28) 



Writing 



Te"^ •/o dx'[w(x')-2XA t (x')f e -^ f* dx'[w(x')-2\A r (x')] 2 = j-^^f* dx' [w{x')-\{A, (x')+A r (x'))] 2 & - f f* dx'[A, (x')-A r (x')] 2 ^ 

(29) 

where the subscripts l,r mean that the operators are respectively to the left and right of V!i{x), we can perform the 
integral over w in ([28]) and, with help of (f20|) . we obtain 



(30) 



= TrV^{x)Te-^^ dx '^ A ^- AR ^ 2 \<j>){(l>\ =TYPg(x)p(t). 
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In (|30|) . Tr is the trace operation and the subscripts L,R mean that the operators are respectively to the left and 
right of =p(0). 

One can proceed further by taking the time derivative of (30). Employing V£(x) = expiH p fP^(x, 0) exp — iH p t 
and (fT7|) . we have 

J^(x) = -iTr{P»(x),H p ]p(t)-± J dxTr[A(x,t),[A(x,t),r£(x)]p(t) 

= -iTr[V$(x),H p )p(t) -\J dzdz'er^-^ 2 Trtf(z, t)£(z,t), [£ f (z', t)&, t), V^x)]p{t). (31) 

The double commutator is easily evaluated. When p, = i, the first commutator is ~ £i(x)£(x), so the double 
commutator vanishes. When p = 0, the first commutator with V p (x) vanishes, so what is left is the double commutator 
with the kinetic energy density T®°(x). This commutator is 

[^Hz,t)^z,t),[^^,t)^,t),^^Hx)-^(x]] =--e t (x)C(x)V 2 5(z-x)-V 2 ^(z'-x). 

Zm m 

Putting this into (|3"Tj) yields the result 
dt 

3A 



p£(a;) = -iTr[P$(x),H p ]p(t)-8' a ^ J dzd Z '5(z-x)6(z' -x)V z -V z ,e-^ {z - z ' }2 TT^(x)ax) P (t) 

-zTr[^(x),i/ p ]p(i) + ^°^ s^^WpW (32) 



Thus, the ensemble-average energy density increase is strictly local, proportional to the expectation value of the 
particle number density. 

We note, when (|32[) is integrated over x, one obtains the well-known result that the particle energy increases linearly 
with time: 



d— 3\N 

dt p Ama 2 



where N is the number of particles. 



VIII. ENSEMBLE-AVERAGE ENERGY-MOMENTUM DENSITY FOR W-FIELD 

Because of the non-local nature of A(x) given in (fT7|) . the local-appearing CQC interaction (|2~Tj) is in fact non-local. 
Thus, as for the particle tensors, the H^-field stress tensor is no longer conserved. A calculation of the expectation 
value of the VF-field stress tensor (equal to the CSL ensemble average) begins by using the expression (11a), 

T^{x) = J dx 1 dx 2 ^[W(x 1 )U(x 2 )-W(x 2 )U(x 1 )]G^(x-x 1 ,x-x 2 ) (33a) 

If 1 

G^ v {x — x\, x — x 2 ) = / dkidk 2 5(kl — k 2 ) 7 s6(s/c°)9(sfc2) sinffci • (x — x±) — k 2 ■ (x — x 2 )] 

(2tt) J ^ i 

\KK + KK + (l/2)r ? ^(fc 1 - k 2 ) 2 ]. (33b) 

The expectation value of (33a) is 

T^(x) = S (y,t\T% v (x,0)\V,t) S = / dxidx2(4>\T f Dwe-^^ dx ' [w{x,) - 2XA ^ )]2 



1 w(x t ) .". ) - w(x 2 ) " ) e-^K dx '^ x '^ 2XA ^'^\(t))G^{x- x u x- x 2 ) 



2 



iSw(x 2 iSw(xi 



= Z! /"* dx X dx 2 {(j)\T J Dwe ~^ ft dx'Mx')-\(Adxl+A r (x'))] 2 e -± J* dx'iAtix^-AAx')] 2 

■[w{x x )A r {x 2 ) - w{x 2 )A r {x 1 )]\<p)G^{x -x x ,x- x 2 ) (34) 
where we have used (|29j) in the second step. Upon integration over w, the result is 

T^{x) = ^ j dx x dx 2 {<j>\Te-i H ^'lMx')-A r (x')f [Al{xi)Ar{x2 ) _ A^A^MG^ (x -x u x- x 2 ). (35) 
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Eq. (|35|) is complicated by the time ordering which embeds Ai , A r within the argument of the exponential. However, 
a simplifying approximation can be made if the particle energies are much larger than a -1 «2eV. When the integrals 
in ([35| over xi, X2 are performed, they act on the gaussians in Ai(xi), A r (xi) and on the argument of the sine in 
G^", giving 

xidx 2 e lkl ' (x ~ Xl) ~ tk2 ' (x ~ X2) e _ ^ (z_Xl)2 e _ ^ (z ' _X2)2 ~ e ''ki-(z-x 1 )- l k 2 -(z'-x 2 ) e -4[ k ?+ k l] . 

That is, k 2 < a~ 2 , approximately. Therefore, we make the approximation, in the expression (33b) for G^ u , 

8{k\ - kl) = 8{kf k 2 - kf + k 2 ) « 5{kf - kf). 
With this approximation, the form factor G in (fT"3f and G 0l/ become 

G(x - sci, 1 - x 2 ) « -<y(x - xi)<5(x - x 2 )e(fi - £2), (36a) 
(^"(s-H.x-sa) « [^ , -^]^(x-x 1 )*(x-x 2 )5(ti-t 2 ). (36b) 
If we put (36b) into (33a), there results 

T° v (x)* ydti[w(x s t)An(x,t)-n(x,t)^-w(x,t)]. (37) 

This (time-independent) expression for the energy-momentum density is what we would have obtained if we had 
looked at the the W-field energy-momentum m (10b), and just removed the integral over x. 
If we put (36b) into ([35]) there results 

Tl U (x) « ^ r^^lTe^/o 1 ^^^)-^^)] 2 ^^)^^^)]^) (38a) 



2 



dt 1 '&[^(x,ti),[A(x 1 t 1 ) J J^]]p(t 1 ). 



(38b) 



In (|38a|) . the exponential's integral's upper limit has been changed from t to ii since the time-ordered integral from 
t\ to t vanishes. The subscripts I, r have been dropped from the commutator since they are evaluated at the same 
time, and are no longer buried within the time ordering of the exponential's argument. We note that, when (|38bj) is 
integrated over x to obtain the W^-field energy-momentum, and the time derivative is taken, the result is the negative 
of the x-integral of (f3"Tj) so the total energy-momentum of particles and H^-field is conserved (the interaction term 
makes no contribution to the ensemble-average energy). 

As in the previous section, one can proceed further by evaluating the double commutator in (|38b[) . obtaining 



TZ (,') = r° 27T3/2 X ma7 I dz(z-^e-M°-») 2 jyt^i^hM^hMh). (39) 



From Q39|) we see that the ensemble-average energy density is negative. It is also non-local, acquiring an increment 
at x, during (t\,t\ + dt\), proportional to the particle number density at time t\ in a volume a 3 around x. 
We also note that the energy in the W-field, according to (|39|) . is 

3XNt 



P°(t) = J dxSfOc) 



4ma 2 ' 

with rate of change equal to the negative of the particle energy rate of change given in the equation following 



IX. PARTICLE CREATION MODEL 



From the last equation in the last section, according to CSL, over the age of the universe, a proton (mass m) gains 
energy and the W^-field loses energy w 10~ 16 mc 2 . This is too small to have any cosmological significance. 

However, one may speculate that collapse could play a significant role in the creation of the universeQ- If the 
universe obeys quantum theory in its initial stages, it is likely that the Hamiltonian which governs it causes its state 
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vector to evolve into a superposition of possible universes. The choice of which universe is actually ours could be due 
to a collapse mechanism. Suppose that CSL provides this mechanism. A significant amount of VF-field energy density 
could be produced, whose gravitational influence could thereafter have a role to play in the further evolution of the 
universe, for example, in its expansion or in galactic formation. 

A simple model of particle creation to illustrate the generation of such M^-field energy has been presented Here, 
we shall reconsider that model, and calculate for it the particle and W-field energy densities. The particle energy 
density is 

H p (x) = m?(x)£(x) + g(^(x) + £(x)], (40) 

so the particle Hamiltonian is H p = J dxH p (x). This describes a "displaced" harmonic oscillator at each x, with 
displacement ~ <?(x). From an initial no-particle state |0) p , without collapse, the ensemble average of the particle 
number density Af p (x) — ^ (x)£ y (x) just oscillates. As we shall see, with collapse, it steadily grows as the oscillation 
dies out. 

We first consider the ensemble-average particle energy density. Eg. (|31[) holds for any Heisenberg particle operator 
density, not just V p (x), so 

-|a7 p (x) = —iTr[f\f p (x), H p ]p(t) — ^ /dxTrL4(x,i), [A(*,t),M p {x)]p{t) 

= -iTr[N p (x), H p ]p{t) = -ig(x)\j?(x) - £(x)], (41) 



^(x) = -iTr£(x), II p ]p(t) - ± I ,/xTr[.l(x./).[.l(x./U(,->]/>(/) 



im + — 



t(x) - i S (x). (42) 



The solution of (gSJ) is 



J( X )= z!gW ri 
U ' im + (\/2) [ 



a _[ im+ (A/2)]t] 



Thus, (|4T|) and (|40|) give respectively 



d_ 



(43) 



(44) 



= ^rw ^ w[1 ~ e " At/2cosmt] ~ r~ At/2sinr 4 (45) 

When there is no collapse (A = 0), we see that Af p (x) just oscillates, and H p (x) remains constant (each oscillator 
at x conserves energy). When A ^ 0, Af p (x) and Ti p (x) asymptotically grow linearly with time, with growth rate 
~ g 2 (x)A. Thus, the permanent creation of particles depends upon the collapse mechanism. 
The W- field energy density can be found from (|38b[) : 

ri 



T w {x) « - J <tt 1 '&[A(x,*i) ) [A(x,ti) J fl r p ]]p(ii) 



'0 

-A 



2(7ra 2 ) 3 / 2 [m 2 + (A/2) : 



J dze-^ {x - z)2 g 2 (z) dii{m[l-e- Atl / 2 cosmi 1 ] - ^ e - A W 2 sin m^}. (46) 



The salient features here are that this energy density is negative, it is generated non-locally (the contribution at x 
comes from a distance as a about x), and it grows linearly with time for t >> A -1 . The total energy is conserved (the 
sum of the spatial integral of the time-derivative of (|4*6")l and the spatial integral of (|43| equals 0) . 



11 



X. W-FIELD ENERGY DENSITY AS A GRAVITATIONAL SOURCE 



We shall consider an example in which the W-ficld energy density is a source of gravitation. The Hamiltonian 
which governs the state vector evolution in the CQC Schrodinger picture is taken to be 



H = J dxT™(x) 



1 - GM 



N v {z) 



M J dxTVp(x) 



GM 



2A J dxA(x)n(x,0) (47) 



In Eq. (|47[) . there are three terms: call them A, B and C respectively. A consists of the W- field energy and the 
gravitational energy of its interaction with particles of mass M. For simplicity, we shall take T®°(x) to be given by 
the approximate expression (|37|) . B consists of the particle mass-energy and the particle gravitational self-interaction 
energy. C is the usual CQC interaction between the W-field and particles, given in (|2"7) . which is responsible for the 
collapse dynamics, where A(x) is given by (|17[) : 



A(x) 



(tto 2 ) 3 / 4 



dze 2 



Mi 



■MJV p (z). 



(48) 



What is not present in (j4"7| is the particle kinetic and potential energy. This is usual when one wishes to display 
collapse behavior without interference by particle dynamics. It is a good approximation when the collapse time is 
shorter than the time scale over which the states evolve appreciably. Also omitted is the W-field gravitational self- 
interaction energy and the gravitational contribution of the particle kinetic and potential energy-densities, which are 
presumed to be negligibly small. 



In the CQC interaction picture, 



from which follows 



\9,t) =e lt{A+B) e- UH \^,0), 



Jt(A+B)(_ 



(-iC)e-"( A+s )|*,t) 
\-tC)e-* tA \*,t), 



(49) 



since B commutes with both A and C. We now expand the initial particle state in eigenstates \n r ) of the number- 
density operator Af p (z) (Af p (z)\n r ) — n r (z)\n r )), 



|*,0> = |0M = |0>5> r ( 



\n r ) 



(for definiteness, the number density eigenstates have been given discrete labels). Since Af p (z) is the only operator in 
the Hamiltonian which acts on the particle state vector, the states \n r ) do not alter during the state vector evolution. 
Thus the operator Af p (z) in A and C in (|4"9"|) may be replaced by the c-number n r (z). With that replacement, A(x) 
given by (|4"8")) will be denoted ^4 r (x). 

To evaluate (|49|) . we need to calculate exp(ityl)n(x, 0) exp(— it A) where 



A 



dx'-[W(x.', <')n(x', t') - n(x', t')w(*.', t')) 



= / dx'[-ii( x \t')w( x ',t')}[i + Mx% 



1 - GM / dz, 



n r (z) 



It follows that 



e l *- 4 n(x, o) e - ltA = n(x, o) + t[i + ^ r (x)]ri(x, o) + . 

Thus, the solution of (gU) is 

|*,t) = ^ Cr (o) e - 42A / ,i * / /o'«« /n (*'.*'[i+*r(x')]Mr(*')| )|, 

r 

= ^C r (0)e^ i2A * dX / O t[1 + '' W1 drn(x,r)A r (x)[l + r (x)] 



n(x,t[i + r ( x )]). 



(50) 



(51) 



(52) 
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In (|52|) we have made a change of variables r = t'[l + <f> r {x')], x = x', so dt'dx' = drdx[l + </> r (x)] _1 . 

We may now proceed, as in (|22[) . to go to the \w) basis and obtain an expansion of the state vector whose terms 
correspond to CSL state vectors: 

|*,t) = ]Tc r (0)K) f Cw | w ) e -^A/ d x/ o t l 1 + ^Wl^,.(x)[l + 0.(xO]- 1 ^(x,r) e - 3 L/ d ^ 2 (,') 

= ]Tc r (0)K) / D ^|^ e -^/^7^d t '[^( x ' >t ')-^(x')[i+^(x')]- 1 e(t[i+0 r (x')]-t')e(*')] 2 

= 5> r (G)|n r ) | I3w | w>/We -i/'ix'/o t[1+ *'' (5t ' ,] ^h( x '-*')-^(x')[i+^(x')]- 1 ] 2 (53) 



where 

rO , ,2/ v ,/\ 



/(*) = 6^ / dX ' /-co ^ 2 (x',t') e -4X / rfx ' + <"(*'>*')_ 

Eq. (|53p is identical in form to Eq. (|22p except for two outstanding distinctions. 

The first is that, in (|22|) , is determined by w(x.',t') values between the hypersurfaces if = and if = t for all 

x'. However, in (|53|) . the relevant u>(x', i') values lie between the hypersurfaces t' = and i' = i[l + r (x')] for all a;'. 
Otherwise, the collapse behavior is the same. The CSL states 

t ) w = J2 c r (0)\n r )e-^ I dx ' * h(x' 1 *')-^(x')[i+<Mx')]-] 2 

T" 

can be plucked from since, according to it, if two w(x', if) are identical between the above specified hypersurfaces 
but differ in the infinitesimal slice hypervolume between t' = (t[l + </> r (x')], (t + dt)[l + 0,.(x')]), the associated states 
are orthogonal forever after. 

Moreover, as is well known, if t is the time read by a clock far from gravitational sources, the gravitational time 
dilation effect is such that, where there is a weak static gravitational potential, a local clock at x' reads time 



if = <V-.9oo(x') = ty/1 + 20(x') « t[l + 0(x')]. 

That is, instead of collapse taking place on evolving t hypersurfaces, the Hamiltonian (|47|) describes collapse as taking 
place on evolving local proper time hypersurfaces. 

The second distinction is that Eq. ([22|) describes collapse toward eigenstates of A(x), but Eq. (|53|) describes collapse 
toward eigenstates of A(x)[l + 0(x)] _1 . Actually, that has no essential effect on the collapse behavior, since A and 4> 
both depend only on Af p (z), and so collapse in both cases is toward the eigenstates \n r ). 

However, there is a subtle difference. The off-diagonal elements of the density matrix have different decay rates. In 
the first case, the decay of the r — s element is proportional to J dx(A r (x) — A s (x)) 2 , while in the second case it is 
proportional to 



dx(A r (x)[l + ^(x)]- 1 - A s (x)[l + <£ r (x)] 



-1\2 



So, the question arises as to how to understand this or, more broadly, how to interpret the appearance of A r (x')[l + 
(MX')]" 1 in Eq. (US]). 

Here is one point of view. Suppose that the gravitational potential of a point mass is not the Newtonian value, but 
rather is that due to the mass smeared by a gaussian over the scale a. So far, there is no experimental restriction on 
this possibility, for a sa 10 -5 cm|9]. Define 

so the gravitational potential is 

0(x) ee -Gm [ dz^t. 

J |x-z| 

According to ([48]) . 

A(x) = (4™ 2 ) 3 / 4 ^ p (x). 
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Now, it is not unreasonable to replace the particle's mass-energy density ~ -A(x) by the particle's total energy density, 
mass-plus gravitational, in the term C in the Hamiltonian, i.e., 

i4(x)-i4(x) = A(x)[l + 0(x)]. 

With the Hamiltonian (|47f replaced by 

H = J dxT°°(x)[l + <Hx)] + M J dxA/" p (x)[l + 0(x)] + 2A j dxA(x)[l + 0(x)]n(x,O), (54) 
the argument goes through as before, with the result 

hM> = $>(0)K> / ^| W )/We^^ x '^ 1+ ' r(5t ' ,Idt ' [M(x '*' ) -^ (x ' )]2 . (55) 



That is, we have put the total energy density -A(x) into the Hamiltonian's interaction term C, but the collapse 
dynamics is determined by A(x), the mass-energy density, as is usual in CSL. 

APPENDIX A: OBTAINING THE ENERGY-MOMENTUM-STRESS DENSITY TENSOR 

Here, a method of obtaining the stress tensor (fTTj) is presented. The idea is to think of W(x) as composed of free 
quantum scalar fields of all possible masses (bradyonic and tachyonic) and energies (positive and negative). Although 
W(x) does not have a Lagrangian, each such quantum field does have a Lagrangian and thereby a stress tensor, and 
when these stress tensors are added up, they yield a stress tensor for W{x). 

Define the quantum fields 

4™,rA*) = (^3/2 / dkS ( k2 + rm 2 )Q(sk a )[b(k)e ik - x + b\k)e~ ik - x ] (Al) 

In (|A1|) . m is the mass of the field, r = 1 or —1 if the field is respectively bradyonic or tachyonic, s = 1 or —1 if the 
energy is respectively positive or negative. One readily finds the commutator 

[4>m,r,s{x),4>m',r',s'{x')) = 6(m 2 - m' 2 )S r r'S ss ' 

• J dke(k 2 +rm 2)[ e ik-(x-xO-i S o; m , r (k)(t-tO + e -ik.(x-x')+i^ m , r .(k)(t-t')] ( A2 ) 



where uj m . r (k) = Vk 2 + rm 2 . For fixed m — m', r = r',s = s', this would be the usual commutation relation, except 
for the factor 8(m 2 — m 2 ) = 6(0). If heuristically we think of 6(0) = 1/dm 2 , we may think of v dm 2 (j) m ^ :S (x) as a 
usual quantum field. Then, we may construct the usual stress tensor for these usual quantum fields, and add them 
up, getting a combined stress tensor for all possible four-momenta, 

T£ v (x) = J2 ^ 2 -{a^ m , r , s (x)a^ m ,, s ( a; ) + 9> m , r , s (x)a^ m , r , s (x) 

rs o 

- jf^[dx<l)m,r, s (x)d x (j) m ^ s (x) + ^(4> m , r ,s(x)d\d x (f> m , r , s (x) + d\d X (j) m ^ s (x)(j> m ^ s (x))] | (A3) 
In (|A3|) . the usual mass-squared term rm 2 (j) 2 n r s (x) has been rewritten utilizing the dynamical equation 

rm 2 <f) m . r ^(x) = d\d x (j) m , r , s (x). 

The factor s in (|A3|) gives the negative energy fields their negative energy. Note that relativistic invariance has been 
broken by treating the stress tensor for the positive and negative energy tachyon fields separately. This is because 
relativistic tachyon quantum field theory requires exchange of tachyon creation and annihilation operators when the 
tachyon energy changes signQ, but no such exchange takes place for b(k), b^(k). 

The four-divergence of this stress tensor vanishes when the dynamical equation is utilized: 

roc 

W = J2 L. ^ 2 4{[^^™,r., S (x)][9A5Vm,r, S ^)] + [9A9Vm,r, S (x)][^0 m , r , s (x)] 

-[d»d X d X 4> m , r , s (x)]4> m .rAx) ~ ^rAxWd X d X <j) m ^ s (x)]} = 0. (A4) 
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That this stress tensor satisfactorily represents the W^-field energy-momentum is shown by calculating the four- 
momentum which follows from (|A3[) . utilizing (|A1[) : 



E 

rs 

E 



dm 2 s 



dkdk'S(k 2 + rm 2 )S(k' 2 + rm 2 )@(sk )Q(sk' )S(k - k')26 t (/c)6(fc) 



{k°k' v + S° v (l/2)[k ■ k' - k' 2 }} 
1 



dfce(s/c )e[r(|fc°| - |k 

|k| 




-|k| 




k°k"2b\k)b(k) 

dk k U b\k)b{k) : 



dkk v tf{k)b(k), 



(A5a) 
(A5b) 

(A5c) 



which is identical to JHJ). 

The term ~ b(k)b(k') (also its Hcrmitian conjugate) is not present in (A5a), for the following reasons. The delta 
function 5(k + k') and step functions make it ~ 6(k, k°)b(— k, k°). When v ^ 0, this is symmetric under change of sign 
and the curly bracketed factor k which multiplies it is antisymmetric so the integral over k vanishes. When v = 0, the 
curly bracketed factor which multiplies it is {— k°k'° + (l/2)[— k ■ k' — k' 2 ]}, which vanishes for k = — k' and fc° = k'°. 

Also in (A5a), b^(k)b(k) + 6(fc)6 T (fc) has been replaced by 2b^{k)b{k). When v ^ 0, the infinite contribution of 
the commutator vanishes since ~ J dkk vanishes. When v = 0, we disregard the infinite vacuum energy as is usual. 
We also note that (A5a)'s curly bracketed term simply becomes k°k u , since the delta and step functions imply that 
k = k'. 

In going from (A5a) to (A5b), we have performed the integrals over k' and m 2 . The integral over k'° of 8(k 2 — k' 2 ) — > 
5[{k° — k'°)(k° + k' )] produces the factor l/(2|fe°|). The integral over m 2 provides the step function limit on the 
range of fc° in (A5b). 

In going from (A5b) to (A5c), we have used ^ s sk°Q(sk°) — \k°\. 

It now remains to express T£ v (x) in terms of the random fields W(x) and H(x), instead of the fields (f> m ,r.s(x). To 
that end we apply I© to (|A1 1) . obtaining: 



1 



(2tt) 7 /2 

and substitute this into (IA3 



dxi J dk5(k 2 + rm 2 )Q{sk°)[^=W{x 1 ) cosfc ■ (x - Xx) - 2\/~XR{xx) sin k ■ (x - xx)], 



(A6) 



T^(x) 



(2 



~yfE^ dm2s J dx l dx * J dk 1 dk28{kl + rm 2 )5{k 2 + rm 2 )Q{skl)Q{skl) 



{ 



\k^k v 2 - (l/2)? 7 ^(fc 1 • k 2 + kl)]\-=W(x 1 )smk 1 ■ (x — Xx) + 2VAII(a;i) cos kx ■ [x — xx] 

L VA 



^-j=W[x 2 ) sinfc 2 • (x - x 2 ) + 2VXl~l{x 2 ) cosfc 2 • (x - x 2 ) 



-{\j2)rf v k\ — -j=W{xx) cos kx • (x — x\ sin kx - (x — xx) 

L\/A 



—=W(x2) cos k% ■ (x — x 2 ) — 2VXIl(x 2 ) sin fc 2 • (x — X2) 
V A 



(A7) 



When m 2 is integrated over and r is summed over, the range of integration of kx, k 2 becomes unrestricted. When 
s and the variables of integration kx , fc 2 are changed in sign in (|A7j) , the result is (|A7j) with the overall sign changed 
and with the sign of the sine's reversed. When half of this is added to half of (|A7|) . the coefficients of W{x\)W{x2) 
and II(a;i)n(a; 2 ) vanish. Moreover, the result is unaltered if we add, to half of it, half the same expression with the 
variables of integration kx, k 2 exchanged and xx, X2 exchanged, obtaining 



^yrE s / dx i dx 2 J dk 1 dk 2 8{k 2 1 - kl)Q{skl)Q{skl) 



2(2tt) 7 

Uk£ kl + k^kg - (l/2)^' y (fci + fc 2 ) 2 ] sin[fci • (a; - xx) + k 2 ■ (x - x 2 )][W {xx)U.{x 2 ) + U{xx)W(x 2 )\ 
+ [k»k» + k\k» + (l/2)r 7 ^(fc 1 - k 2 ) 2 ] sin[fci • (x - x ± ) - k 2 ■ (x - x 2 )}[W(x 1 )U(x 2 ) - Il{xx)W (x 2 )}}.(A8) 
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Now, the first term within the curly brackets in (IA8|1 does not contribute to P£. For it, J dxT®"(x) provides an 
overall factor of <5(k! + k 2 ), so = — k 2 . This, together with the other delta function and step functions gives 
k® = k 2 . Thus, for the first factor in the first term within the curly brackets becomes 

k'lk a 2 + fc°fc 2 ° - (l/2)7 ? 00 (fc 1 + fc 2 ) 2 = 2(fc°) 2 - (1/2)(2A:°) 2 = 

while for P^, this factor becomes 

k-^k<2 ~\~ k\k 2 — ^1(^1 ^1) — 0' 



Therefore, for simplicity we may omit this term, without affecting any of the needed properties of T£ v (x) Moreover, 
one readily sees that 

[W( Xl )Tl{x 2 ) +Tl(x 1 )W(x 2 )} ~ [W + {xx)W + {x 2 )-W-{x 1 )W-{x 2 )\ 

where W+{x) is the positive frequency (annihilation) part of W{x) and W~{x) is its negative frequency (creation) part. 
Thus, by omitting this term, we eliminate the complication of T£ u (x) creating or annihilating pairs of "W-particles." 
This is not the case with the second term we keep, since 

[W{ Xl )IL(x 2 ) - n{x 1 )W{x 2 )\ ~ [W + {x 1 )W-{x 2 ) - W^{x 1 )W + {x 2 )\. 

which maintains the number of "W-particles." 
Therefore, our result is: 

T£ v (x) = J d Xl dx 2 J dk 1 dk 2 5(kl~k 2 2 )Y,se{sk Q 1 )e(sk 2 )sm[k 1 -(x-x 1 )-k 2 -(x~x 2 )} 

[k{k v 2 + k\k^ + (1/2)77^ - fc 2 ) 2 ][^(xi)n(a; 2 ) - n(xi)W(x 2 )]. (A9) 
which is further discussed in section II. 

APPENDIX B: FORM FACTOR 

In this appendix we calculate the form factor which appears in Eq.(llb), which expresses the non-local nature of 
the stress tensor: 

G(x — x\,x — x 2 ) = 1 7 / dk 1 dk 2 8(k\ - /c 2 ) sQ(sk°)Q(sk 2 ) sin[fci • (x - xi) - k 2 • (x - x 2 )]. (Bl) 
(2n)' J , 



To evaluate G, we look at 



7 ( Sl , S2 ) = -i— J dk 1 dk 2 6(kl~ki)e(k a 1 )e^- k ^-\ (B2) 



where Si = [x — x{), and note that 

pOO poo 



2 

g(si,s 2 ) = G(SX,S 2 ) + 



p poo poo 

/ (ikidk 2 (5(fc 2 - fc 2 ) / dkl I dk 2 sin[A:i ■ s% + k 2 ■ s 2 ) 
J J Jo 



so 

1 



G(si,s 2 ) = -\g(si,s 2 ) - g(s 2 ,si)]. (B3) 



To evaluate g, we note that, in (|B2[) , it is 



2iir J n 2 + A 2 



e 



16 



(A is an infinitesimal) which keeps the expression from being relativistically invariant and easy to evaluate. Therefore, 
we take e(fcj) out of the integral: with Tj = s°, we write (|B2I) as 



ff(si,s 2 ) = 







= — / 




2m J 












i 






h 



n 



n 2 + a 2 
n 



i(27r) 7 



dkidkz- 



-Qd/dTi 



n 2 



i{2ir) 8 



2ir 



( ^ a;e i^[fc?-fc2] e i [ fc l- ;i l-' i: 2- s 2] 



where we have written <r = sf — s|, and noted that f2/(f2 2 + A 2 
Thus, we conclude, from ifel) . lIBil) , 



G(si, S2 ) 



(9 2 



(2^) 4 3a 2 



= PI /ft (the principal value). 
Ti6(cr + T 2 ) | T 2 6(-g + T 2 



2 
2 



(B4) 



(B5) 
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